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This paper reports follow-up photometric and spectroscopic observations, includ- 
ing warm Spitzer IRAC photometry of seven white dwarfs from the SDSS with appar- 
ent excess flux in UKIDSS X-band observations. Six of the science targets were selected 
from 16 785 DA star candidates identified either spectroscopically or photometrically 
within SDSS DR7, spatially cross-correlated with HK detections in UKIDSS DR8. 
Thus the selection criteria are completely independent of stellar mass, effective tem- 
perature above 8000 K, and the presence (or absence) of atmospheric metals. The in- 
frared fluxes of one target are compatible with a spatially-unresolved late M or early 
L-type companion, while three stars exhibit excess emissions consistent with warm 
circumstellar dust. These latter targets have spectral energy distributions similar to 
known dusty white dwarfs with high fractional infrared luminosities (thus the i^-band 
excesses). Optical spectroscopy reveals the stars with disk- like excesses are polluted 
with heavy elements, denoting the ongoing accretion of circumstellar material. One of 
the disks exhibits a gaseous component - the fourth reported to date - and orbits a 
relatively cool star, indicating the gas is produced via collisions as opposed to subli- 
mation, supporting the picture of a recent event. The resulting statistics yield a lower 
limit of 0.8% for the fraction dust disks at DA-type white dwarfs with cooling ages 
less than 1 Gyr. Two overall results are noteworthy: all stars whose excess infrared 
emission is consistent with dust are metal-rich; and no stars warmer than 25 000 K are 
found to have this type of excess, despite sufficient sensitivity. 

Key words: circumstellar matter — planetary systems — stars: abundances — stars: 
low mass stars — brown dwarfs — stars: evolution — white dwarfs 



1 INTRODUCTION 

Near-infrared photometry of white dwarfs examines a sci- 
entific phase space that is difficult or impossible to probe 
for main-sequence stars. The first to recognize this potential 
was |Probsd(|l98lh . who realized the compact nature of white 
dwarfs allowed a straightforward photometric detection of 
very cool stellar and substellar companions. Brown dwarfs, 
which often require sophisticated detection techniques when 
searched for at main-sequence stars, are up to 10 times larger 
than a typical white dwarf and can be readily detected in 
spatially-unre solved observations as excess near-infrared fiux 
l|Probstlll983l '). Thus, white dwarfs are an excellent tool for 
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studies of the low-mass stellar and s u bstellar mass function 
via companions (.Farihi et ahl 20051 : IZuckerman fc BecklinI 
1 19921 : iProbst fc O'ConnelllligS^ '). Photometry with Spitzer 
has extended this potential to include closely orbiting brown 
dwarfs an d massive planets although no candidate s are yet 
apparent (|Farihi et al.ll2008bl : [Mullallv et al.ll2007l ). 

Another way in which near-infrared photometry of 
white dwarfs is advantageous is that circumstellar dust or- 
biting within the Roche limit of the stellar remnant be- 
comes heated suffi ciently to emit in this wavelength range 
(|Kilic et al.l |2006|) . This potential is precluded for main- 
sequence stars because the analogous spatial region, where 
km-size or larger solid bodies are tidally destroyed, does not 
extend significantly above their surfaces (Davidsson 199!^). 
Furthermore, any material generated in that narrow region 
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Table 1. Multi- Wavelength Photometry and Derived Parameters for SDSS White Dwarfs 



SDSS Prefix 


0959 


1159"^ 


1221 


1247 


1320 


1506 


1557 


Spectral Type 


DAZ 


DQp 


DAZ 


DA 


DA 


DA 


DAZ 


Spectroscopic Data 


WHT 


SDSS 


WHT 


SDSS 


SDSS 


WHT 


SDSS+WHT 


a (J2000; hms) 


09 59 04.69 


11 59 33.10 


12 21 50.81 


12 47 40.93 


13 20 44.68 


15 06 26.18 


15 57 20.77 


<5 (J2000; ° ' ") 


-02 00 47.6 


+13 00 31.6 


+12 45 13.3 


+10 35 56.1 


+00 18 54.9 


+06 38 45.9 


+09 16 24.6 


FUV (0.15 Mm; ABmag) 




> 20 


19.55 ±0.06 


18.23 + 0.02 


16.52 + 0.03 


18.40 + 0.02 


18.03 + 0.02 


NUV (0.23^111; ABmag) 




19.62 ±0.09 


18.78 + 0.02 


18.33 + 0.01 


16.76 + 0.01 


17.31 + 0.01 


18.16 + 0.01 


u (0.36 /im; ABmag) 


18.52 ±0.02 


18.22 ±0.03 


18.54 ±0.02 


18.63 + 0.02 


17.08 + 0.02 


16.96 + 0.01 


18.50 + 0.02 


g (0.47 /im; ABmag) 


18.09 ±0.03 


18.14 ±0.02 


18.19 ±0.02 


18.50 + 0.01 


17.13 + 0.01 


16.57 + 0.01 


18.45 + 0.01 


r (0.62 /im; ABmag) 


18.36 ±0.02 


17.75 ±0.01 


18.40 ±0.02 


18.82 + 0.02 


17.47 + 0.02 


16.63 + 0.01 


18.81 + 0.01 


i (0.75 /im; ABmag) 


18.54 ±0.02 


17.70 ±0.02 


18.56 ±0.02 


19.13 + 0.02 


17.72 + 0.02 


16.79 + 0.01 


19.13 + 0.01 


2 (0.89 /im; ABmag) 


18.88 ±0.06 


17.82 ±0.02 


18.86 ±0.05 


19.33 + 0.06 


18.05 + 0.03 


16.93 + 0.01 


19.34 + 0.04 


J (1.25 /xm; mag) 


18.32 ±0.06 


17.40 ± 0.03 


18.43 ±0.07 


18.95 + 0.08 


17.62 + 0.04 


16.41 + 0.01 


18.82 + 0.06 


ri (1.64 /^m; mag) 


18.17 ± 0.U9 


17.32 ± 0.U6 


18.39 ± 0.10 


18.90 ± 0.11 


1 T n'-\ 1 n rto 

17.71 ± 0.U8 


16.37 ± 0.02 


19.03 ± 0.14 


K (2.20 /im; mag) 


17.62 ±0.11 


17.08 ±0.08 


18.01 ±0.15 


18.47 ±0.20 


17.53 ±0.11 


16.35 ±0.03 


18.35 ±0.15 


IRACl (3.55 /im; ^tjy) 


80 ±4 


42 ±2 


40 ±5 


25 ±2 


35 ±2 


75 ±4 


35 ±2 


IRAC2 (4.49 /im; ^tjy) 


76 ±4 


28 ±2 


40 ±6 


17 ±2 


25 ±2 


50 ±3 


40 ±2 


Teff (K) 


13280 ± 20 


9500 ± 500 


12250 ± 20 


18540 ± 150 


20220 ± 180 


10670 ± 10 


22810 ± 40 


log g (cms~^) 


8.06 ± 0.03 


8.0 


8.20 ±0.03 


7.86 ± 0.03 


8.37 ±0.04 


8.21 ±0.03 


7.54 ±0.01 


M (M©) 


0.64 ±0.02 


0.59 


0.73 ±0.02 


0.54 ±0.02 


0.86 ±0.02 


0.73 ±0.02 


0.42 ±0.01 


d (pc) 


203 ±4 


110 ±00 


180 ±4 


396 ±9 


150 ±4 


71 ± 2 


566 ± 2 


IR Excess'' 


Dust 


Atm 


Dust 


dM/L 


Bkgd 


None 


Dust 



Note. Ult raviolet, optical, and near-infrared photometry are from GALEX l lMartin et aI.ll2005^ . SDSS l lAbazaiian et al.ll2009l) . and 
UKIDSS l|Lawrence et al.ll2007l 'l respectively. For stars with multiple spectroscopic or photometric datasets, table values are the 
weighted average of available measurements. The SDSS photometry is given in PSF magnitudes; the only catalog entries appropriate 
for point sources. Stellar parameters are derived by fitting the Balmer line profiles as described in §3. 

"This star was selected independently from the six primary science targets (see §3.2). 

'See §3 for the various infrared excess descriptions. 

for 7444, and photometry (only) for 9341, DA star candi- 
dates with g < 19 ABmag. Cross-correlation with UKIDSS 
DR8 then resulted in 1884 objects in common with both 
H- and A'-band detections; these were all fitted with white 
dwarf atmospheric models to probe for A'-band photometric 
excess. Of these 1884 cross-correlated sources, 147 objects 
were identified as having an infrared excess, including 12 ob- 
jects (7 spectroscopic and 5 photometric only) in which the 
excess was most pronounced in the /i'-band, and compatible 
with either a brown dwarf later than L7 or a dust disk; the 
UKIDSS data alone could not typically distinguish between 
those two possibilities. 

Six disk candidate white dwarfs found in this manner 
were selected as targets for Spitzer Cycle 6 observations, 
and are listed in Table [1] Two of thes e targets (1320 an d 
1557) were independently identified bv lSteele et al.l |20l3). 
They p erformed a cross-correla tion of the white dwarf cat - 
alogs of lEisenstein et al.l (|2006l ') and lMcCook fc SionI (|l999l 'l 
with UKIDSS DR8 to identify stars with excess emission 
via optical and infrared colors, as well as spectral fitting. 
A seventh t arget, the DQ peculia r white dwarf LP 494-12 
(SDSS 1159; lEisenstein et ahllioO^ I was selected on t he ba- 
sis of an independent measurement of isT-band excess (jFarihil 
l2009h . corroborated by its UKIDSS photometry. 



2.1 Optical Spectroscopy 

Of the six primary science targets, SDSS spectra were avail- 
able within SDSS DR7 for 1247, 1320, and 1557. In addition, 
long slit spectroscopy was obtained for 0959, 1221, 1506, and 



above the star will rapidly dissipate due to radiation pres- 
sure and drag forces. Thus, while asteroids and comets com- 
monly pass sufficiently close to the Sun, and presumably 
other main-sequence stars that host planetary systems, the 
tidal disruption and subsequent accretion of this debris can 
only be witnessed at w hite dwarfs (and possibly neutron 
stars; IWang et al.ll2006D . 

This paper reports warm Spitzer IRAC measurements 
and optical spectroscopy of six white dwarf candidates se- 
lected for apparent Tf-band excess fluxes based on the 
ground-based photometric surveys the Sloan Digital Sky 
Survey (SDSS) and the UKIRT Infrared Deep Sky Survey 
(UKIDSS). Target selection criteria and survey data are 
given in §2, along with a description of the Spitzer and 
ground-based observations. The spectroscopic and photo- 
metric data analysis is presented in §3, with a summary 
of the overall results in §4. 



2 TARGET SELECTION AND OBSERVATIONS 

The six primary s cience targets were selected as de- 
scribed in detail by ICirven et al.l (|201ll ). Briefly summa- 
rizing, DA white dwarf candidates were selected both pho- 
tometrically and spectroscopically from within SDSS DR7 
l|Abazaiian et al]|2009il . a nd spatially cross-cor related with 
sources in the UKIDSS (|Lawrence et all 120071 ') DR4 (and 
later DR8). A uqriz co lor selection was implemented based 
on the lEisenstein et al. (2006) SDSS DR4 white dwarf cat- 
alog, which was then used to select and download spectra 
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Table 2. Independent JHK Photometry 



4000 4500 5000 

Wavelength [A] 

Figure 1. Normalized SDSS and WHT spectra of the six DA 
white dwarfs in the study are shown in gray, along with the best 
fit models plotted in black. The SDSS spectra stop at 3800 A. 
The WHT spectra of 0959, 1221, and 1557 each contain a narrow 
absorption line near 4481 A, that is due to Mgll. 



Star 


J 


H 






(mag) 


(mag) 


(mag) 


1221" 


18.49 ± 0.05 


18.21 ± 0.05 


18.05 ±0.05 


1221'' 


18.50 lb 0.05 


18.27 ±0.05 


18.22 ±0.10 


1320" 


17.64 ±0.05 


17.44 ±0.05 


17.49 ±0.05 


1320'' 


17.60 ±0.05 


17.61 ± 0.07 


17.58 ±0.11 


1557 


19.05 ±0.05 


18.92 ±0.05 


18.56 ±0.05 



"Aperture photometry 
'PSF-fitting photometry 



APALL within IRAF. Finally, the instrumental response was 
removed by dividing each observation by the response func- 
tion, calculated by dividing the associated standard by its 
corresponding flux table. 



1557 during two observing runs in 2010 April and May us- 
ing using ISIS, the dual beam spectrograph mounted on the 
4.2 m William Herschel Telescope (WHT). The primary aim 
of these observations was to corroborate the photometrically 
selected DA candidates (0959, 1221, 1506) and to improve 
on the SDSS spectroscopy of 1557. For this purpose, the 
R600 grating was used in the blue arm with a l" slit and a 
2 (spectral) by 3 (spatial) pixel binning. This setup covered 
the wavelength range 3650 — 5120 A, i.e., the entire Balmer 
series except Ha, with an average dispersion of 0.88 A per 
binned pixel in the blue arm, and 7691 - 9184 A with 0.99 A 
per binned pixel in the red arm. The red spectra had rela- 
tively low signal-to-noise ratio (S/N). 

The blue spectra were de-biased and flat-fielded using 
the STARLlNpQ packages KAPPA and FIGA RO and then opti- 
mally extracted using the PAMEL^^code ()Marshlll989[ ). The 
extracted spectra were wavelength calibrated using CuNe 
and CuAr arc lamp exposures, and finally fiux calibrated 
using observations of appropriate standard stars obtained 
with the same instrumental setup. The photometric classifi- 
cation of 1221, 1506, and 1557 as hydrogen dominated white 
dwarfs was confirmed by the ISIS spectroscopy. Normalized 
WHT and SDSS spectra of all Spitzer DA targets are shown 
in Figure [T] 

Additional observations of 09 59 and 1320 were obtai ned 
with the VLT using X-Shooter l|D'Odorico et all |200d ) in 
service mode in 2010 June- July and 2011 February-March. 
The raw frames were reduced using the X-Shooter pipeline 
version 1.3.7 within GASGANC0. The standard recipes were 
used with default settings to extract and wavelength cali- 
brate each spectrum. The final extraction of the science and 
spectrophotometric standard spectra was carried out using 



^ Developed and maintained by the Joint Astronomy Centre and 
available from http://starlink.jach.hawaii.edu/starlink 
^ Developed and maintained by T. R. Marsh and available from 
jhttp: / /www. War wick, ac.uk/go/ trmarsh 
http:// www.eso.org/sci/software/gasgano 



2.2 Additional Near-Infrared Photometry 

Independent near-infrared photometry for 1221, 1320, and 
1557 was obtained on the 23 March 2011 with the WHT 
using the Long-Sl it Intermediate Resolu tion Infrared Spec- 
trograph (LIRIS: iManchado et al.1ll998h . Images taken in 
a 9-point dither pattern were obtained in the J-, and 
ifs-band filters (MKO system; R. Karjalainen 2011, private 
communication) with typical total exposure tim es of 270 s in 
clear conditions. Three standard star (ARNICA: lHunt et al.l 
119981 ) fields were observed in a similar manner for photo- 
metric zero-point calibration. The data were reduced in the 
standard manner, by subtracting a median sky from each 
image in the dithered stack, flat-fielding (using sky fiats), 
then averaging and recombining frames. 

LIRIS suffers from what is known as a detector reset 
anomaly, which appears in certain frames as a discontin- 
uous jump (in dark current) between the upper and the 
lower two quadrants. To remove this unwanted signal, af- 
ter flat-fielding and sky subtraction, the detector rows were 
collapsed into a median column (with real sources rejected) , 
and subsequently subtracted from the entire two dimen- 
sional image. The resulting fully reduced frames exhibit 
smooth backgrounds, free of the anomalous gradient. 

Aperture photometry of standard stars was performed 
using r = 3'.' 75 aperture radii and sky annuli of 5" — 7'.' 5 
in size. For the relatively faint science targets, smaller pho- 
tometric apertures were employed with corrections derived 
from several brighter stars within the same image field and 
filter. For both 1221 and 1320, point-spread function (PSF) 
fitting (i.e., DAOPHOt) was used in addition to photometry 
with small apertures. All data taken in the J^s-band filter 
were fiux-calibrated using the ARNICA if-band standard 
star photometry, and the error introduced by this should be 
significantly smaller than the 5% absolute calibration un- 
certainty. The independent JHKs photometry is listed in 
Tabled 
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SDSS 0959 
T = 13280 K 



SDSS 1231 
T = 12250 K 



0.010 r: 



0.001 



1.0 10.0 



0.001 I 



1.0 10.0 
X (^tm) 



1.000 F 



0.100 r 



0.010 - 



0.001 L 
0.1 



SDSS 1247 

T = 1B540 K 



1.0 10.0 



1.000 



0.100 



0.010 



0.001 I 

100.0 0.1 



SDSS 1320 

T = 20320 K 



1.0 10.0 



1.000 F 



0.100 - 



0.001 U 
0.1 



A (/xm) 



SDSS 1506 
T = 10670 K 



0.001 I 

0.1 



SDSS 1557 
T = 228 10 K 



1.0 10.0 
A (/xm) 



Figure 2. Ultraviolet through infrared spectral energy distributions of the DA-type, infrared excess candidates. Stellar atmosphere 
models are plotted as dotted lines, using parameters derived from model fits to hydrogen Balmer lines in the SDSS or WHT spectra, 
and matched to the observed g-band fluxes. Table [T] photometry is shown as data points with error bars. 



2.3 Spitzer IRAC Observations 

Near-infrared imaging observations of the white dwarf tar- 
gets were obtained with the warm Spitzer Space Telescope 
(| Werner et alj 120041 ) during Cycle 6 u sing the Infrared Ar- 
ray Camera (IRAC; iFazio et al] 120041 ') at 3.6 and 4.5 /im. 
The total integration time in each channel was 1200 s, where 
the observations consisted of 40 frames taken in the cycling 
(medium) dither pattern with 30 s individual exposures. All 



images were analyzed and fluxes measured as in lFarihi et alJ 
(|2010bl ) using 0.6 pixel mosaics created using MOPEX. In 
cases where the flux of a neighboring source was a potential 
contaminant of the white dwarf photometry (i.e., 1221 and 
1247, see §3.1), steps were taken to minimize or remove any 
such external contributions, including small aperture pho- 
tometry and PSF fitting with daophot and apex. 

Of all IRAC targets, only 1506 is detected in the 
WISE Preliminary Data Release. Reliable data exists only 
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Table 3. Relative Number Abundances in Disk-Polluted Stars 



SDSS 1557 (WHT) ^ 



SDSS 0959 (X-shooter) 



SDSS 0959 (WHT) 



0.5 -Tj '/v'^'^'lfv^VyV' ^^yV-vJ^^y^V ^- 

" I SDSS 1281 (WHT) 



3900 3920 3940 3960 4440 4460 4480 4500 4520 

Wavelength [A] 

Figure 3. Black lines are model fits to the Call 3934 A and Mgll 
4481 A line regions in the normalized ISIS and X-Shooter spectra, 
plotted in gray, of all three white dwarfs found to have circum- 
stellar dust. The stars are displayed top-to-bottom by decreasing 
eflfective temperature: 22 810 K, 13 280 K, 12 250 K. 



at 3.4 /im, where the flux (76 ± 8/iJy) is in agreement with 
the IRAC photometry, albeit with larger errors. Generally, 
the science targets are too fai nt for WISE, whe re confusion 
can be a serious concern (see lMelis et al.ll201ll '). and hence 
the choice of warm Spitzer. 



3 RESULTS AND ANALYSIS 

The SDSS and WHT spectra of the six DA white dwarfs were 
fitted using mode l atmosp here spectra computed with the 
code described by Kocstcr ( 2010'), and fo llowing the method 
described by Jlebassa-Mauscrgag (|2007l ). In brief, the grid 
of model spectra was fitted to the normalized Balmer line 
profiles, leading to a 'hot' and 'cold' solution of roughly equal 
absorption line equivalent width. This degeneracy is broken 
by fitting the models also to the slope of the flux-calibrated 
spectra. The best-fit models are overplotted on the data in 
Figures [T] and (2] and the parameters are reported in Table 

m 

All targets stars but one reveal excess emission at IRAC 
wavelengths when compared to the atmospheric models, de- 
rived from their optical spectroscopy. However, there are 
only four genuine cases where the excess radiation is asso- 
ciated with an additional component physically associated 
with the system: three dust disks and one (likely) low-mass 
stellar companion. All individual stars are discussed below 
in some detail. 



3.1 New Disk- Polluted White Dwarfs 

Three stars - 0959, 1221, 1557 - display infrared fluxes 
consistent with dust disks orbiting within their respective 
Roche limits for large asteroids. The excess emissions are 
too strong, and the infrared col ors are incorrect for sub - 
stellar or planetary companions (|Farihi et al.ll2009l . l2008al ). 
Specifically, even for the case of a substellar object with 
anomalously red colors, the 4.5 flux requires its size to be 



Object 


log(Mg/H) 


log (Ca/H) 


Mg/Ca 


Sun 


-4.45 


-5.66 


16 


Chondrites 


-0.72 


-1.96 


17 


0959 


-5.2 


-7.0 


59 


1221 


< -5.6 


-7.5 


< 85 


1557 


-4.5 


< -5.7 


> 15 



Note. Measured abundance uncertainties are 0.2 dex. 



several to ten times that of Jupiter, in contrast with observa- 
tions and models of brown dwarfs and planets [Leggett et al.l 
(2010). Firmly corroborating the disk interpretation, each of 
these DA white dwarfs was found to be metal-lined in fol- 
low up WHT and VLT optical spectroscopy, consistent with 
ongoing infall of circumstellar material. Specifically, their 
ISIS spectra contain narrow absorption lines near 3934 A 
and 4481 A, which are due to CallK and Mgll, respectively; 
they are DAZ stars. Figure displays model fits to the the 
detect ed metal absorption features, using tlusty and SYN- 
SPEC iHubenv fc Lanj|l995l : lLanz fc Hubenvlll995l ) in order 
to estimate the photospheric abundances. 

Only Mg and Ca are detected, and Table [3] reports the 
determined abundances and upper limits for these two el- 
ements. The Mg/Ca ratio observed in 0959 is significantly 
larger than that found in chondrites, but similar to that 
found for the highly polluted white dwarf GALEX 1931 
ijVennes et al.ll2011i ). where the relative poverty of Ca, Ti, 
and Al is potentially consistent with the removal of crust 
and mantle in a differentiated, terrestrial-like body during 
the p ost-main sequence evolution of the star (|Melis et al.l 
I2OIII ). The Mg abundances determined from the X-Shooter 
and ISIS spectra of 0959 are consistent to within the errors, 
while the low S/N around the CallK line prevents a mean- 
ingful abundance estimate for that element based on the ISIS 
data. For 1221 and 1557, only upper limits are derived for 
Mg and Ca respectively, partly due to the inferior quality of 
the ISIS spectra compared to the X-Shooter data. However, 
another factor is that the undetected atomic transitions are 
not easily e xcited at the respective effective temperature of 
these stars (jZuckerman et al. 2003). 

For all three stars, the excess infrared e missions ar e fit- 
ted with optically thick, flat disk models (| Jural [200^ ). As 
discussed in previous white dwarf-disk studies, optically thin 
dust models are inferior, pri marily due to the Foynting- 
Robertson timescales involved (|von Hippel et arii2007.) . Not 
only would such dust be removed on the timescale of years, 
but the particles orbit once every few hours, ensuring the 
disk will relax rapidly into a flat conflguration; a geometri- 
cally thin, optically thick disk circumvents these issues. Im- 
portantly, such a disk can harbor sufficient dust mass to ac- 
count for the most highly polluted helium atmosphere stars 
Farihi et ai]|2010al ). whereas an optically thin disk cannot 



Jura et al.ll2009l ) 



The fitted disk parameters are listed in Table [4l and 
the model fluxes are shown together with the observational 
data in Figure [4] For each star, R/d is derived from the at- 
mospheric models using the spectroscopically derived stellar 
parameters and observed photometry: R is determined from 
Tcft and log g while d is specifled by m — M for the same 
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Table 4. Disk Parameters 



Star 


Tcff 


R/d 


Tin 


'-pa 




' out 


cos i 




(K) 


(10-12) 


(K) 


(K) 








0959 


13280 


1.37 


1600 


800 


10 


25 


1.0 


1221 


12250 


1.41 


1400 


800 


11 


23 


0.7 


1557 


22810 


0.73 


1400 


800 


25 


52 


0.5 



'^Owing to the wavelength coverage, the outer disk temperature 
and radius is not as well-constrained by the data as the inner disk 
temperature. See §3.1 for a detailed discussion of uncertainties in 
the disk parameters. 

(e.g., g — Mg, r — Mr). For fixed R/d, there are three free pa- 
rameters for a flat disk; inner disk edge temperature, outer 
disk edge temperature, and inchnation. These models con- 
tain a modest amount of degeneracy in their ability to fit 
the data, even when longer wav elength data, suc h as 8 and 
24 /im photometry are available jjura et al.ll2007l ). However, 
the models provide good, if broad constraints on the geom- 
etry and temperature of the disks. The inner disk edge is 
fairly well constrained by the excess emission at 2.2 and 
3.6 fim, but not uniquely so, while the 4.5 fim fiux can be 
reproduced by a relatively wide temperature range and a 
higher inclination, or a narrower temperature range and a 
lower inclination. This is because the total disk emission is a 
function of its temperature profile and solid angle, and the 
latter is comprised of the inclination and the ring annuli. 

Despite this uncertainty, the infrared excesses cannot be 
satisfactorily fit with inner disk temperatures below 1200 K, 
nor with significant amounts of T < 500 K dust. These new 
discoveries essentially mimic the thermal emission profiles 
observed at 16 dusty white dwarfs with longer wavelength 
data in addition to t he wavelengths covered by warm Spitzer 
l|Farihi et al.ll2010bl '). The excess at 0959 is sufficiently strong 
that it requires a nearly face-on disk, while the disk incli- 
nations for 1221 and 1557 are less constrained. Importantly, 
the excesses are consistent with inner disk temperatures ap- 
proaching or exceeding that which causes rapid sublimation 
of typical dust grains; gas thus produced gives rise to viscous 
drag an d enhances in fiow of disk material onto the stellar 
surface (|Rafikovll201lh . While the outer disk temperatures 
are not more tightly bound, the range of acceptable values 
places all the debris within 1.5 i?©, where large, solid bod- 
ies will be shredded. Overall, the data are consistent with 
circumstellar debris originating in large, star-grazing solid 
bodies that is now being gradually falling onto the stellar 
surface. 



3.1.1 0959 

This white dwarf displays Can emission lines that are the 
hallmark of closely-orbiting metallic gas. These emission fea- 
tures are clearly detected in the X-shooter data, and, upon 
close inspection, are marginally visible in the noisy red ISIS 
spectrum. Both red spectra are shown in Figure O along 
with SDSS 1228, the prototypical white dwarf with emitting 
metal gas and dust (|Brinkworth et cil]|2009l : iGansicke et al.l 
l2006l ). Notably, the Call line profiles for this star are sub- 
stantially narrower and weaker compared to those detected 
in the other three published cases. In 0959 the full width at 
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Figure 4. Infrared excesses fit by circumstellar dust disk models 
whose parameters are listed in Table |4] The plot features are 
the same as in Figure [J] but on a linear scale to emphasize the 
infrared excesses. Dashed lines represent emission from optically 
thick, flat disk models, and a the dashed-dotted lines represent 
the sum of the stellar and circumstellar model fluxes. 



zero intensity and equivalent width are close to 250 km s ^ 
and 3 A, respectively, while the known cases ha ve values in 
the range 1000 - 1400 knts"! a i id 10 - 60 A (|Mehs et all 
I2OIOI : iGSnsicke et~al]|2008l . |2007| . 120061 ). Assuming that the 
structure of the gas disks in 0959 and 1228 are similar, the 
velocity width of the Call lines in the two systems can be 
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Figure 5. The Call triplet region of 0959 in the normalized 
WHT and X-Shooter spectra exhibit clear emission lines from 
a gaseous disk component. For comparison, also plotted is the 
normalized spectrum of the prototype gas disk around SDSS 1228 
llGansicke et al.l2006l) . The narrower emission features at 0959 in- 
dicate a disk that is significantly inclined relative to that seen at 
1228, and consistent with the high fractional infrared luminosity 
measured with IRAC. 



used to estimate the inclination of the disk orbiting 0959. 
For 12 28, two independent estimates yield an inchnation 
of 70° l|Melis et al.ll2010l : iBrinkworth et al.l 120091 ) ■ and this 
yields i ~ 12° for 0959. The independently estimated disk 
inclination from the model fitted to the infrared emission 
is 10°, consistent wit h a near face-on disk of dust and gas. 
iGansicke et al.l (|2006h found that the Ca ll emission lines in 
1228 are optically thick, and hence a low disk inclination 
will naturally result also in a lower equivalent width. 

SDSS 0959 becomes the fourth published case of a 
white dwarf with detectable, gaseous debris in addition 
to solid dust particles in a circumstellar disk ( Melis et al.l 
l2010l : lFarihT et al ] |2010bl : IBrinkworth et al ]l200gl). At Teff = 
13 300 K it is by far the coolest star to host such dual-phase 
debris, demonstrating conclusively that the detected gas in 
dust disks is not related to grain sublimation in the presence 
of higher temperature white dwarfs. In the three well-studied 
cases, the gaseous and solid debri s are essentially spa - 
tially coincident (Melis et al. .2010. : Brinkworth et"ani2009h . 



both spanning a regio n from around 20 to 100 stellar radii 
(jCansicke et al.ll2006l ). Critically, there is a distinct lack of 
emission from gas in the innermost regions where sublima- 
tion is expected, and which may be due to lower disk surface 
density, resulting from e nhanced viscosity and inward drift 
of gas (jMelis et al.ll201ol ). Corroborating this picture, sub- 
limation of dust at 0959 is expected only within 15 stellar 
radii, and hence cannot account for emitting gas out to 100 
radii. 

It is uncertain if the current disk at 0959 and other gas- 
dust disk hosts are recently created or extant disks that have 
been impacted by a additional body or material. Regardless, 
the presence of detectable gas is almost certainly related to 
the destruction of solid material. 



3.1.2 1221 

The IRAC images of this star reveal a point-like source 
neighboring the white dwarf. Both sources were simultane- 
ously modeled and deconvolved with PSF-fitting photome- 
try. The neighbor is faintly but clearly visible in the UKIDSS 
-ff-band image, and is identified as a separate source in that 
survey. The LIRIS J^s-band images place it 1'.'8 distant dis- 
tant at position angle 256°. Together with its IRAC fluxes, 
the JHK photometry of the neighb or yields colors inco nsis- 
tent with a (substellar) companion (|Patten et al.ll200^ '). Al- 
though relatively faint, it appears to be diffuse and extended 
in the ground-based images and is thus a background galaxy. 
Both the UKIDSS and LIRIS photometry of the white dwarf 
reveal strong infrared excess at K, with the latter dataset 
confidently free of photometric contamination by the neigh- 
bor, in agreement with the excess determined from IRAC 
photometry. No Can lines are detected in the low S/N red 
ISIS spectrum. 



3.1.3 1557 

The mass of the white dwarf in 1557 i s significantly below 
the average mass of field white dwarfs fci ebert et al.ll2005l ). 
While the fit to the SDSS spectrum is consistent with a low- 
mass, carbon-oxygen core, the WHT data are more sugges- 
tive of a helium core. The SDSS data are of relatively poor 
quality at S/N ^ 10, while the WHT spectrum has S/N 
~ 30. Table [T] lists the weighted average of these two solu- 
tions, favoring the low mass interpretation. If a low surface 
gravity is corroborated in future observations, this may sig- 
nify an innermost planet that was consumed during the first 
ascent giant phase, ejecting suffic ient envelope mass to pre - 
vent the onset of helium ignition (|Nelemans fc Taurij|l99^ '). 
This picture is consistent with the evidence presented here 
for a remnant planetary system. As in the case of 1221, no 
Can emission is seen in the poor quality ISIS data. 



3.2 Unusual Near-Infrared Colors for 1159 



This DQ peculiar star (LP 494-12, WD 1156-1-132) was se- 
lected on the basis of independent JHK photo metry that 
revealed a 0.3 mag excess at K band (jFarihii200 9l) relative to 
model expectations for a 10 000 K , helium atmosphere w hite 
dwarf (jHolberg fc BergeronI lioO^ : iFontaine et al.ll200lf l: its 
UKIDSS photometry independently corroborates this result. 
Due to the C2 absorption features at blue-green wavelengths, 
a pure helium atmospheric model was fitted only to the 
izJH photometric data, as shown in Figure |6l The IRAC 
fluxes reveal only a Ray leigh- Jeans tail to the spectral en- 
ergy distribution, but at a shifted, lower temperature rela- 
tive to the 0.6 — 1.5 /j m continuu m. As found for another 
DQp star (LHS2293; iFarihil |2009l ) , this apparent excess is 
likely a re-distribution of emergent flux due to the highly 
absorbed regions within the atmospheric C2 bands. Thus 
the infrared excess only exists relative to a pure helium at- 
mosphere, and is likely due to a distinct, absorption- induced 
shape in the energy distribution of the star. 
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Figure 6. Spectral energy distribution of the DQ peculiar star 
LP 494-12, which ha-s H - K > 0.25 mag in both UKIDSS and 
independently obtained photometry. Fitted to the izJH photom- 
etry is a 9500 K stellar model for a pure helium atmosphere white 
dwarf (see §3.2). 



3.3 A Likely Ultracool Companion at 1247 

The infrared excess observed at this star is distinct, as the 
3.6 /im flux excess is strong, but the 4.5 flux is consis- 
tent with a Rayleigh- Jeans type slope (see Figure [2]). At first 
glance this may seem puzzling, but most M dw arfs approach 
this s pectral behavior at I RAC wavelengths (jPatten et alJ 
l2006l : ICushing et al.l [ioosi ). Neither the UKIDSS nor the 
IRAC images indicate an additional source (i.e., a spatially 
detectable companion), implying any secondary would have 
a projected separation smaller than around 70 AU, and es- 
sentially ruling out a background object as the cause of the 
observed emission. While the infrared excess is solid, and 
possibly apparent at J band, the S/N of the UKIDSS pho- 
tometry is below 10 at H and K, making it hard to esti- 
mate a companion spectral type based on these data. As- 
suming the observed excess emanates from a companion at 
the correct distance, its 2.2, 3.6, and 4.5 /im absolute bright- 
nesses are 11.2, 10.0, and 9.8 mag. Comparing these values 
to the empirical absolute magnitudes of ultracool dwarfs at 
these w avelengths yie lds an approximate spect ral type of 
LO ± 1 dLeggett et al.|[2 010: Patten ct al. 20o d: IVrba et al.1 
l2004l : lDahn et all 12003 : iKirkpatrick fc McCarthvlll994l ) for 
a low mass companion. Superior JHK photometry would 
tighten up this estimate, and possibly indicate a new L-type 
companion to a w hite dwarf, which are rare ({Steele et al.1 
I2OIII : iFarihi et all 2005). Radial velocity monitoring of the 



white dwarf has the potential to detect variability if the sys- 
tem is a close binary, whereas high angular resolution imag- 
ing may be required if the system has a moderately wide 
separation. 



3.4 A Likely Background Galaxy at 1320 

The IRAC images of this star have full widths at half max- 
imum that are inconsistent with a single point source; they 
are enlarged by about 20% in both channels and there is 
modest elongation in the 3.6 /im image roughly along P.A. 
45°. Nothing is seen in the UKIDSS images, but a second 



source appears faintly in the LIRIS H and Kg-hand images 
at a position consistent with the IRAC data. The image of 
this additional source is diffuse and extended and almost 
certainly a background galaxy. Moreover, there is an off- 
set between the image centroids of the science target in the 
LIRIS and IRAC images, indicating that the IRAC source 
(and hence flux) is coincident with the additional source im- 
aged with LIRIS, and not the white dwarf. The flux from 
this likely galaxy falls within I'.'O of the white dwarf, and 
contaminates all available infrared photometry. Even with 
PSF-fltting photometry on the LIRIS images, the resulting 
colors are slightly too red for a 20 200 K white dwarf. Civen 
the overlapping proximity of this source to the white dwarf, 
the infrared excess is probably due entirely to this back- 
ground object. 



3.5 No Near- Infrared Excess for 1506 

Somewhat ironically, the only star without a corroborated 
infrared excess via IRAC photometry is the brightest source 
in the sample, and thus the star with the highest S/N 
UKIDSS photometry. However, this star was selected for the 
Spitzer program when the evaluation criteria were still tied 
to UKIDSS DR4, and has since dropped out of the selec- 
tion based on the r elease of DR8 and an improved analysis 
(|Girven et aLlbOHl ). The model for 1506 shown in Figure H 
shows that an excess is neither observed with Spitzer, nor 
in if -band. While this white dwarf was observed with ISIS, 
conflrming its nature as a DA star, it is hence dropped from 
further analysis. 



4 DISCUSSION AND CONCLUSIONS 

4.1 An Unbiased Dust Disk Frequency 

What makes this study so far unique is that previous 
IRAC programs have either targeted known metal-rich white 
dwarfs, or near-infrared bright white dwarfs. While the flrst 
approach has been highly successful and acco unts for the 
bulk of white dwarf disk dis coveries to date (|Farihi et al.l 
l2010bl . [2OO9I : lJura et al.ll2007l '). two drawbacks exist. First, 
it requires 8 m class telescope time using high- resolution 
spectroscopy to efficiently identify (weak) metal lines i n 
white dwarfs l|Koester et al.l l2005l : IZuckerman et al"] |2003| ) , 
and second, only about 20% of metal-contaminated white 

J '■ " 'n 

dwarfs have disks fFarihi et al. 2009). The second method, 
employed by the iMuUallv et al. (2007.) survey of more than 
130 white dwarfs with Ks < 15 mag, was unbiased by the 
presence or absence of atm ospheric metals, but o nly two 
stars with dust were found (Ivon Hippel et al]|2007h . and is 
hence only 1.5% efficient by number of targets. Here, five 
candidate DA white dwarfs with if-band excess have been 
selected with a success rate of 60% for circumstellar dust 
and 80% for confirmed infrared excesses physically associ- 
ated with the system. One of five targets has an infrared 
excess probably due to a background galaxy. 

Based on the c andidate selection c riteria and numbers 
given in Table 6 of lCirven erall l|201ll) . of 1884 candidate 
and confirmed DA white dwarfs with the necessary UKIDSS 
detections, there are 12 stars whose data are consistent with 
disk-like, if-band excesses. Several years of Spitzer studies 
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have shown that only 50% of dusty w hite dwarfs have such 
is'-band excesses iFarihi et al.ll2010bl ). implying that a de- 
cent estimate of the disk fraction among DA white dwarfs 
is at least 1.2%. However, as seen in this study, in some 
cases the A'-band data can be the result of neighboring back- 
ground objects, as well as real companions such as low mass 
stars or brown dwarfs. Thus a more realistic estimate is 3/5 
of this number, implying at least 0.8% of DA (and presum- 
ably all) white dwarfs have dust disks (an d atmospheric met - 
als), consistent with previous estimates (jFarihi et al.ll2009h . 
These results only apply to stars with cooling ages less than 
1 Gyr, as this corresponds to the approximate cut off in ef- 
fective temperature fSOOOK: iFontaine et al.l l200lh for the 
DA color selection. However, Spitzer studies have shown 
a distinct lack of dust disks at older and cooler metal- 
contaminated stars; only G166-58 has an (anomalous) in- 
frared excess and a cooling age beyond 1 Gyr (Farihi et al.i 
l2008bl ). 



4.2 A Lack of A'-Band Emitting Dust Disks at 
Warm to Hot White Dwarfs 

Interestingly, the DA selection criteria is quite sensitive to 
warm to hot, hydrogen-rich white dwarfs, as the Balmer 
decrem ent gives these stars unique colors in that wavelength 
region (jGirven et al.ll201ll ). However, there are no candi- 
dates for dust among thousands of stars(!). Either dust disks 
at white dwarfs with Teff > 25 000 K exist yet rarely ex- 
hibit A-band excess or dust at such stars is itself rare. The 
latter possibility is somewhat contradictory based on the 
paradigm of planetary systems that are dyna mically rejuve- 
nated at the onset of the post-main sequence ( Bonsor e t al.l 
I2OIII : [ Bonsor fc Wvattll20ld : iDebes fc Sigurds son 2002). If 
this apparent dearth of warmer systems at A-band is real, 
it would likely require a physical mechanism to preclude 
such emission (such as rapid removal of warm grains, or 
long timcscales for dust disk spreading toward the star). 

Farihi (2011,) plot the fractional infrared luminosity of 
all known white dwarfs with dust circa mid-2010 as a func- 
tion of cooling age, revealing a potential trend of increasing 
dust emission towards cooler stars. This could be a natural 
consequence of the potential for cooler white dwarfs to host 
wider disks: as a star cools and its luminosity decreases, dust 
grains may persist closer to the stellar surface prior to sub- 
limation. Thus, any disk extending inward from the Roche 
limit has the potential to be wider at cooler white dwarfs. 
Conversely, at increasingly warmer stars, the inner dust disk 
edge will be physically further from the star as the radius 
at which grains are rapidly sublimated increases (see Table 
|4]), eventually exceeding th e Roche limit for large solid bod- 
ies l|von Hippel et al] 120071 ). While it is uncertain at which 
precise effective temperature an evolving dust disk might be 
completely sublimated, the results of this work clearly show 
that disk emission analogous to G29-38 is rare among warm 
and hot white dwarfs. More sensitive surveys with Spitzer 
and WISE are needed to address this matter more compre- 
hensively. 



ACKNOWLEDGMENTS 

The authors thank the referee Sandy Leggett for feedback 
which improved the quality and clarity of the manuscript. 
This work is based in part on observations made with the 
Spitzer Space Telescope, which is operated by the Jet Propul- 
sion Laboratory, California Institute of Technology under a 
contract with NASA. Some data presented herein are part 
of the Sloan Digital Sky Survey, which is managed by the 
Astrophysical Research Consortium for the Participating In- 
stitutions (http://www.sdss.org/), and the UKIRT Infrared 
Deep Sky Survey. 



REFERENCES 

Abazajian, K. N., et al. 2009, ApJS, 182, 543 
Bonsor, A., Wyatt, M. 2010, MNRAS, 409, 1631 
Bonsor, A., Mustill, A. J. Wyatt, M. 2011, MNRAS, 414, 
930 

Brinkworth, C. S., Gansicke, B. T., Marsh, T. R., Hoard, 

D. W., Tappert, C. 2009, ApJ, 696, 1402 
Gushing, M. C, Rayner, J. T., Vacca, W. D. 2005, ApJ, 

623, 1115 

Davidsson, B. J. R. 1999, Icarus, 142, 525 
Dahn, C. C, et al. 2002, AJ, 124, 1170 
Debes J. H., Sigurdsson S. 2002, ApJ, 572, 556 
D'Odorico, S. et al. 2006, SPIE, 6269, 626933 
Eisenstein D. J., et al. 2006, AJ, 132, 676 
Farihi, J. 2009, MNRAS, 398, 2091 

Farihi, J. 2011, in White Dwarf Atmospheres and Circum- 
stellar Environments, ed. D. W. Hoard (Berlin: Wiley- 
VCH) 

Farihi, J., Beckhn, E. E., Zuckerman, B. 2005, ApJS, 161, 
394 

Farihi, J., Becklin, E. E., Zuckerman, B. 2008a, ApJ, 681, 
1470 

Farihi J., Barstow M. A., Redfield S., Dufour P., Hambly 

N. C. 2010a, MNRAS, 404, 2123 
Farihi, J., Jura, M., Lee, J. E., Zuckerman, B. 2010b, ApJ, 

714, 1386 

Farihi, J., Jura, M., Zuckerman, B. 2009, ApJ, 694, 805 
Farihi, J., Zuckerman, B., Beckhn, E. E. 2008b, ApJ, 674, 
431 

Fazio, G. G., et al. 2004, ApJS, 154, 10 
Fontaine, G., Brassard, P., Bergeron, P. 2001, PASP, 113, 
409 

Gansicke, B. T., Koester, D., Marsh, T. R., Rebassa- 
Mansergas, A., Southworth J. 2008, MNRAS, 391, L103 

Gansicke, B. T., Marsh, T. R., Southworth, J. 2007, MN- 
RAS, 380, L35 

Gansicke, B. T., Marsh, T. R., Southworth, J., Rebassa- 

Mansergas, A. 2006, Science, 314, 1908 
Girven, J., Gansicke, B. T., Steeghs, D., Koester, D. 2011, 

MNRAS, 417, 1210 
Holberg, J. B., Bergeron, P. 2006, AJ,132, 1221 
Hubeny, I., Lanz, T. 1995, ApJ, 439, 905 
Hunt, L. K., Mannucci, F., Testi, L., Migliorini, S., Stanga, 

R. M., Baffa, C, Lisi, F., Vanzi, L. 1998, AJ, 115, 2594 
Jura, M. 2003, ApJ, 584, L91 

Jura, M., Farihi, J., Zuckerman, B. 2007, ApJ, 663, 1285 
Jura, M., Farihi, J., Zuckerman, B. 2009, AJ, 137, 3191 



© 0000 RAS, MNRAS 000, 000-000 



10 J. Farihi et al. 



Kilic, M., von Hippel, T., Leggett, S. K., Winget, D. E. 

2006, ApJ, 646, 474 
Kirkpatrick, J. D., McCarthy, D. W. 1994, AJ, 107, 333 
Koester D. 2010, In Mcmoric dclla Societa Astronomica 

Italiana, 81, 921 
Koester, D., Rollenhagen, K., Napiwotzki, R., Voss, B., 

Christlieb, N., Homeier, D., Reimers, D. 2005a, A&A, 432, 

1025 

Lanz, T., Hubeny, I. 1995, ApJ, 439, 875 

Lawrence, A., et al. 2007, MNRAS, 379, 1599 

Leggett, S. K., et al. 2010, ApJ, 710, 1627 

Liebert J., Bergeron R, Holberg J. B. 2005, ApJS, 156, 47 

Manchado, A., et al. 1998, SPIE, 3354, 448 

Marsh, T. R. 1989, PASP, 101, 1032 

Martin, D. C, et al. 2005, ApJ, 619, LI 

McCook G. P., Sion E. M. 1999, ApJS, 121, 1 

Mehs, C, Jura, M., Albert, L., Klein, B., Zuckerman, B. 

2010, ApJ, 722, 1078 
Melis, C, Farihi, J., Dufuour, P., Zuckerman, B., Bur- 

gasser, A. J., Bergeron, P., Bochanski, J., Simcoe, R. 2011, 

ApJ, 732, 90 

Mullally, F., Kilic, M., Reach, W. T., Kuchner, M., von 
Hippel, T., Burrows, A., Winget, D. E. 2007, ApJS, 171, 
206 

Nelemans, G.,Tauris, T. M. 1998, A&A, 335, L85 
Patten, B. M., et al. 2006, ApJ, 651, 502 
Probst, R. 1981, Ph.D. Thesis, University of Virginia 
Probst, R. 1983, ApJS, 53, 335 

Probst, R. G., O'Connell, R. W. 1982, ApJ, 252, L69 
Rafikov, R. R. 2011, MNRAS, 416, L55 

Rebassa-Mansergas, A., Gansicke, B. T., Rodn'guez-Gil, P., 
Schreiber, M. R., Koester, D. 2007, MNRAS, 382, 1377 

Steele, P. R., Burleigh, M. R., Dobbie, P. D., Jameson, R. 
F., Barstow, M. A., Satterthwaite, R. P. 2011, MNRAS, 
416, 2768 

Venues, S., Kawka, A., Nemeth, P. 2011, MNRAS, 413, 

2545 

von Hippel T., Kuchner M. J., Kilic M., Mullally F., Reach 

W. T. 2007, ApJ, 662, 544 
Vrba, P., et al. 2004, AJ, 127, 2948 

Wang, Z., Chakrabarty, D., Kaplan, D. L. 2006, Nature, 
440, 772 

Werner, M. W., et al. 2004, ApJS, 154, 1 
Zuckerman, B., Becklin, E. E. 1992, ApJ, 386, 260 
Zuckerman, B., Koester, D., Reid, I. N., Hiinsch, M. 2003, 
ApJ, 596, 477 



© 0000 RAS, MNRAS 000, 000-000 



